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Abstract. Stationary electric transport in semiconductor nanostructures is stud- 
ied by the method of nonequilibrium Green functions. In the case of sequential 
tunneling the results are compared with density matrix theory, providing almost 
identical results. Nevertheless, the method of Green functions is easier to handle 
due to the availability of an absolute energy scale. It is demonstrated, that the 
transport in complicated structures, like quantum cascade lasers, can be described 
in reasonable agreement with experiment. 

1 Introduction 

Transport through semiconductor nanostructures Jl],^]|| is dominated by quan- 
tum effects. A corresponding quantum transport theory can be based on 
density matrices or nonequilibrium Green functions (see, e.g., Refs. and 
references cited therein). In this article, I want to show that nonequilibrium 
Green functions provide a strong tool to handle this situation for stationary 
transport, even for complicated structures like the quantum cascade laser. By 
a direct comparison with the density matrix method the differences between 
both methods are examined for the test case of sequential tunneling between 
neighboring quantum wells. 

The article is organized as follows: After a general formulation of the prob- 
lem (Section |^) it will be shown in Section ^ how the method of nonequilib- 
rium Green functions can be applied in the stationary state. A simple exam- 
ple is discussed in detail in Section ^, and a direct comparison with density 
matrix theory is made in Section |^. Finally, the full power of the Green func- 
tion approach is demonstrated in a simulation of quantum cascade lasers in 
Section ^|. 

2 General aspects of quantum transport 

The starting point for a quantum kinetic description is the Hamilton operator 
in second quantisation (using basis states labeled by a) 



H — H a + U + H scatt 



(1) 
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where (in Heisenberg representation) 



H = J2 E <* a t( t ) a «( t ) 



(2) 



a 



is diagonal in the basis a), 



(3) 



describes nondiagonal parts of the Hamiltonian as well as the presence of 
electric fields, and H sca ±t refers to interactions with phonons, random impu- 
rity potentials (which are treated within impurity averaging), or interactions 
between the particles. The final goal is to calculate various observables such 
as the occupation of the state a 



between the respective states, which can be obtained from the equation of 
continuity for the occupations. Here (...) denotes the quantum mechanical 
expectation value with the (nonequilibrium) distribution. Besides these one- 
particle density matrices (ag(t)a a (t)) , higher order density matrices describe 
correlation effects and response functions. 

Two formalisms exist to treat the quantum problem [i.e. to find approx- 
imations in order to to obtain solutions for Eqs. (U]|)]: Within the method 
of density matrices, the temporal evolution of these quantities (where all 
operators are taken at the same time) is studied directly. This method was 
extremely successful in the study of electron kinetics on short time scales, 
see, e.g., Ref. J5| for details. On the other hand, Green functions depend on 
two different times. These have been used for stationary transport and for 
electron kinetics, see, e.g., Refs. |,| for details. 

3 Method of Nonequilibrium Green Functions 

The key quantities in the theory of nonequilibrium Green functions are the 
correlation function (or 'lesser' Green function) 



which describes the occupation of the states (for equal times and indices), 
together with the respective correlations both in time and state index, as well 
as the retarded Green function 



f a (t) = (4 (i)a„(t)> 



(4) 



or transition rates 




(5) 



G< ua2 {h,t 2 )=i{ai 2 {t 2 )a ai {t 1 )) 



(6) 



'ret 
act ,ct2 



(ti,t 2 ) 



iO(h - t 2 )(a ai (h)ai 2 (t 2 ) + al 2 (t 2 )a ai (h)) 



(7) 
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which describes the response of the system at time t\ in state a\ after an 
excitation at time t 2 in state a 2 - 

If the external potential is static and transients resulting from initial 
conditions have disappeared, the system is typically in a stationary statejj] and 
all functions depend only on the time difference t\ — t 2 . Then it is convenient 
to work in Fourier space defined by 



F, 



!1)a2 (£) = ~ / dte iEt / h F aua2 (t + t 2 ,t 2 ). (8) 



This defines the energy E which is not the level energy E a of a certain state, 
but a new parameter, setting an absolute scale to compare energies of different 
states. Then the following equations determine the Green functions 



(9) 







Gii ,« 2 (E) = GtAE)E< , (E)Gp: a2 (E) (10) 

where the advanced Green function is just given by G^(E) = ^G^piE)^ . 

These equations contain the self energies, which are typically functionals of 
the Green functions, and depend on the approximation chosen. To provide a 
glimpse of their structure, two examples are given here: 

In the simple Born approximation, the retarded self energy for impurity 
scattering with scattering matrix elements V aii p is given by: 



S ata 2 {E) = ^2(V aij /3V 0ta2 ) imp E _ E ^ + , Q+ (11) 

where (V aij) 3V(3 ia2 )i m p denotes the averaging over all impurity configurations. 
A comparison with Fermi's golden rule gives 9{^°^(i? Q )} = —h/2T a , which 
relates the imaginary part of S ret to the lifetime T a of the state. 

In the self-consistent Born approximation, the lesser self-energy for phonon 
scattering reads 

P,l PlA 

x [n B (f^i(p))G^ 2 (E - Rw,(p)) ( 12 ) 
+ [n B (M(p)) + 1] Gf u/32 (E + Rw,(p))l , 



1 This is not neccessarily the case as self-sustained oscillations or even chaotic 
behavior frequently occur in semiconductor systems at high electric fields M. 
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Fig. 1. Sketch of the 
test structure consist- 
ing of two neighbor- 
ing quantum wells. 
A fc-conserving tran- 
sition is indicated 



where T (p, I) is the matrix element of the electron-phonon interaction. 
p and I denote the wave vector and the mode of the phonon spectrum with 
corresponding frequency wj(p). n B (E) = (exp(E/k B T) — l) -1 is the Bose 
distribution, describing the occupation of phonon modes in thermal equilib- 
rium. Taking into account that G < describes the occupation of states, one 
can identify the first term as phonon absorption from states with energy 
E — Tiuui (p) , and the second term as the combination of stimulated and spon- 
taneous phonon emission from states with energy E + Huji (p) . This demon- 
strates that E< % a2 (E) is just the in-scattering rate at energy E. Furthermore 
it becomes clear, that E is the relevant energy scale for the occupation, and 
indeed in thermal equilibrium one finds 

G< a (E)=iA a (E)n F (E-v) (13) 

with the spectral function A a (E) = — 23{G„ ct Q ,(£')} and the electro-chemical 
potential /i. It is crucial to note, that the energy E enters the Fermi function 
np(E — u) instead of the level energy E a . 

Eqs. (|9|,|l0|) have to be solved self-consistently with the functionals for the 
self-energies, which is a formidable task. Finally, the physical observables of 
interest can be calculated via the density matrices 

(4(t)a a (t)) = -i J ^G< p (E). (14) 

Such calculations have been performed for the double barrier resonant tun- 
neling diode || and superlattice structures |lO|,[ll]]. In Section ^, simulation 
results for a quantum cascade laser structure are presented. In order to study 
the general concept, an analytic solution for a simple example will be given 
in the subsequent section. 

4 A simple example 

Let us consider sequential tunneling between two neighboring quantum wells, 
labeled by a and b, see Fig. [|. In both wells only the lowest bound state (with 
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energy Eq) plays a role and the in-plane behavior is determined by plane 
waves with quasimomentum k = [k x , k y ) and the dispersion E k — h 2 k 2 /2m. 
Such an experiment has been performed by applying separate contacts to the 
individual layers |l^,^3|. Reflecting the experimental situation, the following 
assumptions are used in the theory: 

• The translational invariance of the structure allows only for /c-conserving 
transitions a — > b with a tunneling matrix-element M. 

• Due to weak coupling between the layers, only a small current is present, 
so that both layers are in thermal equilibrium with the respective contact. 

• The electron densities in both layers are equal, so that the bias U deter- 
mines both the difference in electro-chemical potentials /i a — /i& = elf, 
and the shift of the energy level b by — ell. 

• Impurity scattering is the dominant scattering process, which is realistic 
for the case of vanishing temperature, T — > 0, considered here. Correla- 
tions in the impurity potential between both wells are neglected. 

First note, that the simple application of Fermi's golden rule forbids tran- 
sitions between the wells (see the arrow in Fig. for U =/= 0, due to the 
energy conserving (5-function. Therefore a quantum transport calculation is 
indispensable. From the assumptions mentioned above we obtain the Hamil- 
tonian 



H = Y^^o + E k)a{a k + {Eq + E k - eU)b{b k + Mb{a k + Ma\b k 

k 



(15) 



fe.fe' 



where the matrix elements V kk i for impurity scattering are subject to impu- 
rity averaging during the calculation. 

Now we apply the Green function approach discussed in the last sec- 
tion. The general state index a is replaced by the two indices a/b,k. Ne- 
glecting correlations in the impurity potential between both wells implies 
(Vkk'^k' ft) imp ~ $n,u- Therefore the self-energies are diagonal in the well 
index and Eq. @ takes the form 

Gl e :(k;E) Glf(k;E) 

E-E -E k -E^(k;E) -M 

-M E-E Q -E k + eU - S l b f (k; E) 

(16) 

which gives G rct (fc; E) = G TCt ^(k; E) + AfG rct W(fc; E) + 0(M 2 ) with 



(17) 
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Fig. 2. Transitions 
in the Green-function 
treatment at fixed 
energy E. Transi- 
tions are allowed 
if A a (k,E) and 
At,(k,E) overlap in 
the energy window 
between the electro- 
chemical potentials. 
The inset depicts 
the current-voltage 
characteristic 



and 



G ret(1) (fc;£) 



G^ {0 \k;E)Glf°\k;E)" 




(18) 



Then Eq. @ yields 

G< a (k;E) =G^ \k;E)MG<^Hk ] E) + G^°\k;E)MG^(k;E) 



0{M 2 



(19) 



The zeroth order terms refer to the uncoupled layers in equilibrium. Therefore 
Eq. (O) applies and with Eqs. (J5|Jl4|) one obains the total current M 



(20) 



ja — >b 



2(f or Spin)^-V / ^A a (k; E)A b (k; E) 

k 

x [n F (E - p.) -n F {E- fx b )] + 0(M 3 ) . 



This result shows that the transitions occur at fixed energy E, see Fig. |[ 
The spectral functions A a (k, E) and A b (k, E) give the weight of the levels a 
and b at the energy E, respectively. For each energy, there is a contribution 
to the current if there is a finite weight of both states and difference in 
occupation [np(E — fi a ) — n F (E — fJ-b)] - This demontrates the crucial role of 
E as an energy scale. Following Sec. 3.3 of Ref. |ll[] this expression can be 
simplified for 17* {k; E) = -iP/2: 



ra — >b 
1 GFT 



A 



e\M\ 2 PQ 



2r 



{eUf + r 2 



J dE [n F (E - n a ) - n F (E - /x 6 )] (21) 



where po = m/(nh 2 ) is the two-dimensional density of states including spin 
and A is the sample area. Eq. ( pl| ) was used in the interpretation of the 
experiments P^|L3| where good agreement was found. 
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5 Treatment by density matrices 

For comparison, the example of the preceding section will now be treated 
within the density matrix approach. Here the key quantities are the oc- 
cupations f a {k) — (a\.a k ) and fb(k) = {b\b k ) as well as the polarisation 
P(k) = (b\.a k ), which provides the transition rates according to Eq. (||). The 
key idea is to derive the temporal evolution of these quantities as a set of 
differential equations and to solve these after some approximations. Thus 

~ t P(k,t)=([H,bl(t)a k (t)}) 

= - eUP(k, t) + M [f a (k, t) - f b (k, t)] (22) 
+ E [ V k>k(bl(t)a k (t)) - V k a k ,(bl(t)a k ,(t)) 

k' 

The density matrices {b\{t)a k > it)) contain the phase of the scattering matrix- 
element Vu k i ■ Thus one evaluates the temporal evolution as 

~Vk\>{b\{t)a k/ {t)) = (E k -eU- E k ,)V k a k , (b{a k ,) 

+M[v k a k ,(aia k ,)-V k a k ,(bib k ,)} ( 23 ) 

+ E V k"kV k % ( b l a ^ - E V ™ V *v> 

k" k" 

After impurity averaging one finds (V k „ k V kk ,) imp = | V k \, \ 2 5 va 8 k >^ k > and we 
have 

~V k a k ,(bl(t)a k ,(t)) ={E k -eU- E k ,)V k a k , (b{a k ,) - \V^\ 2 P(k,t) 
i dt (24) 



M 



V k a k ,(aia k ,)-V k a k ,(b{b k ,) 



In Eq ( p4| ) it is tempting to neglect the last line, as these terms contain both 
M and V. In this case one obtains 

V k %ma,(t)) = Ek _ eU " Ek , +i0+ \V k %'\ 2 m (25) 

in the stationary state, where the term i0 + ensures that correlations vanish 
for t — > — oo. Therefore 

Y,Kk>(bl(t)a k >(t)) = Z I a et (k,E + E k -eU)P(k) (26) 

k' 

where Eq. (pT|) has been used. In the same way one finds J2 k ' ^fc'fcC'l' (t) a k(t)) = 



£% dv (k,E + E k )P(k). Inserting into Eq. @ gives the polarisation 



M[f a (k)-f h (k)} 

eU + S^(k, E + E k - eU) - Ef Y (k, E + E k ) 



P(k) = _ TT , ^ _T ^ V „ , ^ (27) 
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and with ^ dv (fc, E Q + E k ) - ^(k, E Q + E k - eU) « if the current 



^DMTsimp = 2 ( for Spin)e— (pTT)2 



[fa(k)-fb(k)} 



(28) 



This means that the transitions occur between states with the same k as 
indicated in Fig. [j]. In contrast to the treatment by Fermi's golden rule, life- 
time broadening allows for these transitions even for finite U. This structure 
appears frequently in density matrix calculations, see, e.g., Rcf. (IJ] for se- 
quential tunneling. On first sight ioMT-simp seems to be similar to Iq^^. But 
for the situation considered, ^MTsimp vanishes, because f a {k) = fb(k) as the 
electron density is identical in both wells. This indicates that the result is 
not trustworthy. 

Now we take into account the first term in the last line of Eq (24), which 
indeed is of the same order as the other terms because P(k) is of order M. 
Therefore we have to evaluate 



Tdt 



ni'KWVW) =(Ek - E k ,)V k a k ,(ala k , 



+M 



[V kk ,{bla k ,)-V k \,{a{b kl ) 
■J2v k a „ k V k a k ,{al„a k r) - V k a k ,V k a lk „{ala k „) 



(29) 



Perfoming impurity averaging and neglecting the M terms (yielding expres- 
sions of order M 2 ) gives: 



V k a k ,(a k (tfa k ,(t)) 



-1 



E k - E w + i0+ 
Now the solution of Eq ( p4| ) provides 

1 



V k %\ 2 [fa(k') - /„(*:)] 



(30) 



V k \.(bl(t)a k ,{t)) 



E k -eU- E v + ^ 



+ 



M 
eU 



1 



1 



E k 

ra 1 2 



eU-E k 



E k - E k , + M+ 



(31) 



x\Vk a k ,\ 2 [fa{k')-Uk)]. 



Using S{l/(x + i0 + )} = — tt5(x) and neglecting the corresponding real part 
one obtains 



YVk\'(b\{t)ak>(t)) « K et (k,E + E k - eU)P(k) 

k' 

M 

+ i— ^{ST(k, E + E k - eU)} [n a (E k - eU) - n a (E k ) 



(32) 
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Fig. 3. Transitions in 
the density matrix 
treatment according 
to Eq. @. The cur- 
rent is driven by im- 
purity assisted tran- 
sitions between states 
of different k, so that 
the energy is con- 
served 



where n a (E k ) = f a (k). In the same way one obtains 
£^*< 6 U*K(*)> » Z% dv (k,E + E k )P(k) 

k' 

M A 

+ i—ME^ik, E + E k )} [n b {E k ) - n b {E k + ell)} 
eu 



(33) 



Inserting into Eq. (^) and using ,r; adv / rct ~ ±_%r/2 one obtains after a few 
lines of algebra 



ra — >b 
i DMT 



2r 



h ^ {euf + r 2 



n a (E k ) - n b (E k + eU) n a (E k - eU) - n b (E k ) 



(34) 



This means, that the current is driven by the difference in occupation at 
the same energy, rather than the difference in occupation at the same k. 
Fig. ^ suggests that these transitions should be viewed as scattering-assisted 
transitions. Note, that this interpretation is completely different from the 
result by Eq. (^). The fact, that the transitions occur at fixed energy strongly 
resembles the Green-function formalism, see Fig. ||, and the evaluation of 
Eq. (^f|) and Eq. (HJ) give similar results. The difficulty in obtaining this 
result in density matrix theory seems to be related to the fact that, unlike in 
the Green function formalism, the energy scale E is not available. 



6 Results for quantum cascade laser structures 

Electrically driven unipolar semiconductor lasers can be made by a proper 
design of nanostructures. The key point is population inversion between sub- 
bands, which can be achieved by a rather complicated sequence of quantum 
wells. In order to amplify this effect, several periods of such sequences are 
grown on top of each other, resulting in a quantum cascade laser [ |l5| oper- 
ating in the 10 |im range. 
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Fig. 4. Conduction 
band of the quantum 
cascade laser struc- 
ture from Ref. @. 
The eight lowest 
Wannier functions 
are given for one 
period, where the 
base line depicts 
the respective level 
energy. Note that 
these functions are 
repeated in each 
period 

In Fig. U such a structure (from Ref. |L6|]) is shown. Here one period 
consists of an active region (where the lasing transition occurs) of 3 wells 
(0 < z < 15 nm) and the injector region of 5 wells (for 15nm < z < d). This 
sequence with period d = 45.3 nm is repeated 30 times in the experiment. 

Similar to the treatment in Section^, we use a set of basis states consisting 
of localized states in the growth direction (labeled by v) and a free particle 
behavior with wave vector k perpendicular to the growth direction. For the 
localized states we use Wannier functions as depicted in Fig. [I| This allows 
for a calculation of all matrix elements in Eqs. (|^,|3|) from the sample parame- 
ters following Ref. [jll). For field strengths of the order of F w 50 kV/cm, the 
states in the injector become aligned to the excited state (dashed-dotted line) 
in the active region of the next period. This design provides a strong enhance- 
ment of population of this level associated with with population inversion in 
the active region. 

The self-energies are evaluated within the self-consistent Born approxima- 
tion treating scattering with impurities, optic phonons, and acoustic phonons 
similar to Ref. [jnj. Eqs. (p[jIo|) were solved self-consistently with the function- 
al for the self-energies reaching an accuracy of approimately 1%. Details will 
be given elsewhere [jlTj . For practical purposes the momentum dependence of 
the matrix elements had to be neglected. Thus, these quantum transport cal- 
culations can be considered as complementary to semiclassical Monte-Carlo 
simulations |l^,^9|, where such details can be easily taken into account. 

Fig. shows some results from our calculations. The field-current relation 
is in good agreement with the experimental findings of Ref. Jl6| . The electron 
densities of level v are evaluated by 

k 



A 




-40 -20 20 40 
z [nm] 



Transport in nanostructures 11 



1 1 1 1 — 

Level 1 

Level 2 

Level 4 
Level 5 
Level 7 




1 i i i i 

2 4 6 
J [kA/cm 2 ] 



Fig. 5. Voltage drop per pe- 
riod Fd as well as occupation 
n v for different levels versus 
current density. The linestyle 
corresponds to the fat lines 
in Fig. ^j. The occupation of 
the other levels is below 5 x 
10 10 /cm 2 for all currents 



For vanishing current and field, essentially the lowest level is occupied. With 
increasing field, the energy of the levels 4, 5, and 6 is diminished because 
they are located further to the right. Thus, the carriers are transfered to 
these states, which are located in the injector region. At the same time the 
occupation of level 7 increases linearly with the current, which indicates that 
the transport essentially occurs via this level. (The sample was designed pre- 
cisely for this behavior.) At moderate currents, occupation inversion occurs 
between this level 7 and the levels 1, 2, and 3, which are lower in energy. As 
the dipole matrix element between these levels and level 7 is rather large (all 
of them are essentially located in the active region), lasing is likely to occur, 
which is indeed observed in the device for current densities above 7 kA/cm 2 . 

7 Conclusion 

Nonequilibrium Green functions provide a powerful tool to study quantum 
transport in nanostructures under stationary conditions. They allow for a 
quantitative description of nanostructure devices such as superlattices and 
quantum cascade lasers. The power of this method results from the availabil- 
ity of the global energy scale E, the frequency related to the time difference 
in G{t\, t 2 ). Although essentially the same results can be obtained in density 
matrix theory, the structure of the equations appear more complicated there. 

Fruitful cooperation with S.-C. Lee and financial support by DFG within 
SFB296 and FOR394 are gratefully acknowledged. 
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